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SUMMARY 


A systematic study is made of the approximate inviscid theory of 
thin bodies moving at such high supersonic speeds that nonlinearity is 
an essential feature of the equations of flow. The first-order small- 
disturbance equations are derived for three-dimensional motions involv- 
ing shock waves, and estimates are obtained for the order of error 
involved in the approximation. The hypersonic similarity rule of Tsien 
and Hayes, and Hayes' unsteady analogy appear in the course of the devel- 
opment . 

It is shown that the hypersonic theory can be interpreted so that 
it applies also in the range of linearized supersonic flow theory. 

Hence, a single small-disturbance theory, and associated similarity 
rule, apply at all supersonic speeds above the transonic zone. 

Several examples are solved according to the small-disturbance 
theory, and compared with the full solutions when available. These 
include flow past a wedge and cone, and determination of the initial 
gradients at the tip of plane and axially symmetric ogives. For the 
ajcially symmetric ogive it is shown that further terms can be found 
only by using Lighthill's technique of rendering solutions uniformly 
valid, and thus the initial curvature of the pressure distribution is 
calculated. It is concluded that on a body of revolution described by 
a power series, the pressure distribution and shock wave are also given 
by power series. 

A brief discussion is given of various additional approximations 
from existing theories. 


INTRODUCTION 


Aerodynamic shapes are ordinarily most efficient when they cause 
the least flow disturbance. For this reason, simplified theories based 
upon the assumption of small disturbances due to thin bodies have proved 
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to be of practical value in analyzing incompressible, subsonic, tran- 
sonic, and supersonic flows.^ For flows at Mach nxambers large compared 
with unity, however, the pressure disturbances may no longer be small 
(compared with the static pressure) even for thin shapes, so that in 
this sense it has been said that no small-disturbance theory exists 
(ref. l). However, if viscosity can be neglected, the velocity disturb- 
ances remain small compared with the speed of flight (though not compared 
with the speed of sound), and even the pressure changes are small if com- 
pared with the dynamic pressure. In this sense, therefore, a small- 
disturbance theory exists, and the assumption of such small disturbances 
leads to a useful simplification of the equations for compressible flow 
at arbitrarily high Mach numbers. 

Viscosity and heat conduction must be neglected in order to have a 
small-distiirbance theory. Otherwise, for example, the viscous no-slip 
condition would introduce velocity disturbances equal to the speed of 
flight. In many cases this simplification does not destroy the essential 
features of the flow. In other cases, the inviscid theory may serve as 
a basis for including viscosity and heat conduction. Thus, recent studies 
of the hypersonic boundary layer (refs. 2 and 3)^ which indicate that vis- 
cous effects become essential at extreme. Mach numbers (say, greater than 
15), replace the boundary layer by a fictitious solid surface and then 
utilize inviscid theory of the sort considered here. 

At sufficiently high Mach mmibers, inviscid flow past any given thin 
object requires nonlinear equations for its description. We take this as 
the definition of hypersonic flow: Supersonic flow past a thin body is 

termed hypersonic if the Mach nmber is so great that nonlinearity becomes 
an essential feature. Thus, the definition of hypersonic flow stands on 
an equal footing with the generally accepted meaning of transonic flow 
at the other extreme of the supersonic range; that is, flow at a Mach 
number so close to unity that nonlinearity (of a different sort) is an 
essential feature. These two terms - transonic and hypersonic - are most 
meaningful when defined (as here) only for thin shapes. They then describe 
two quite distinct regimes which are, moreover, separated by a consider- 
able range of "ordinary supersonic" flow in which the transonic and hyper- 
sonic nonlinearities are unimportant, so that linearized theory can 
account for all significant features of the flow. If one attempts to 
extend the terms to thick bodies, the two separate regimes tend to merge, 
so that one must concede that a flow field can be simultaneously tran- 
sonic and hypersonic. 

It should be noted that the term hypersonic has occasionally been 
used in the literature with other meanings than that adopted here. Flows 
are sometimes called hypersonic if the free-stream Mach number is simply 
large compared with unity (say, 10 or or even 3)* This is not strictly 


throughout, "thin" is used to refer to any body whose streamwise 
slope is small, and so applies to slender fusiform objects as well as 
flat shapes such as airfoils. 
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equivalent to the present definition because, in principle, at any given 
Mach number a body can always be chosen so thin that nonlinearity is 
insignificant. However, such extreme thinness does not arise in prac- 
tice, so that the two definitions are equivalent for practical purposes. 
Again, Oswatitsch has defined hypersonic flow as the limiting condition 
for a given body as the free-stream Mach n\miber tends to infinity 
(ref. 4 ). This is a limiting case of the present definition and, indeed, 
Oswatitsch 's similarity rules for thin bodies are simply special cases 
of the more general rxiles . 

Associated with each of the various small-disturbance theories is 
a similarity rule which connects flows at different speeds past affinely 
related shapes. In the case of linearized subsonic and supersonic 
theory, the similarity rule was fully understood only long after the 
small-disturbance theory was in common use. On the other hand, the tran- 
sonic similarity rule was developed concurrently with the small- 
disturbance theory. For hypersonic flows, the simplified theory and 
associated similarity rule were first given by Tsien (ref. ^), but were 
restricted to irrotational flows (and to plane or eixially symmetric 
shapes ) . This is a severe limitation because strong curved shock waves 
and consequent entropy gradients and flow rotation are- essential features 
of nearly all hypersonic flow problems. This restriction was removed by 
Hayes, who indicated in a brief note (ref. 6 ) that Tsien’ s similarity 
rule is valid for rotational flows and for general three-dimensional 
shapes . The rule was further extended to unsteady motion by Hamaker and 
Wong (ref. 7 ). 

As a result of this circuitous development, there exists a gap in 
the hypersonic theory. The similarity rule is known for full three- 
dimensional flows with curved shock waves, but the underlying small- 
disturbance theory has never been written down. (To be sure, however, 
its form is known from the analogy with nonsteady flow in one less 
dimension, which was pointed out by Hayes.) For the special case of 
plane flow, and for Mach numbers which are not arbitrarily large, this 
gap has recently been closed by Goldsworthy (ref. 8 ). The published 
examples of applications of the theory are limited to the few special 
cases which are strictly irrotational (e.g., the wedge and cone) or are 
assumed to be approximately so. 

The present paper undertakes a systematic study of the small- 
disturbance theory for hypersonic flow. First, the small -disturbance 
problem is derived by reduction of the full equations of motion, bound- 
ary conditions, and shock-wave relations. The similarity rule and 
unsteady analogy appear in the course of this development. This portion 
of the paper may be regarded as an elucidation of Hayes’ note, with esti- 
mates obtained for the order of error. Next, it is pointed out that to 
within terms of the order neglected, the hypersonic similarity theory 
can be written in the form of the similarity theory for linearized super- 
sonic flow, so that a single theory and associated similarity rule cover 
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both regimes. Then, a nvunber of special problems are solved according 
to the small -disturbance theory and are compared with the full solutions 
when they exist. 

The symbols' used in the text are defined in Appendix A. 

HYPERSONIC SMALL-DISTURBANCE THEORY 
Basic Assumptions 

Consider a three-dimensional body fixed in a steady uniform stream. 
Viscosity and heat conduction are neglected, which implies that shock 
waves will be approximated by abrupt discontinuities. 

The body is assumed to be thin, in the sense that the streamwise 
slope of its surface is everywhere small compared with unity. The degree 
of thinness will be measured by the small parameter t which may , for 
example, be taken to be the maximum slope of the body,^ or its thickness 
ratio. However, for inclined shapes x must be identified with the 
angle of attack if it is considerably greater than the body thickness. 

The free-stream Mach number M is assumed to be so high that the 
flow is hypersonic. That is, linearized theory is inadequate for pre- 
dicting the essential features of the flow. It is known that linearized 
theory yields an adequate approximation if the maximum body slope is 
small compared with the slope of the free-stream Mach cone, that is, 

if pr«l where P = As this ratio approaches unity, linear- 

ized theory grows increasingly inaccurate. Therefore the flow is hyper- 
sonic if the ratio px is not small compared with xmity. Since x is 
small, this means that p and, therefore, also M will be large in the 
hypersonic range, so that p is nearly equal to M. Thus, the criterion 
for hypersonic small-disturbance flow may be written 


X « 1 ]. 
M » 1 j 


with Mx ~ 1 or » 1 


( 1 ) 


From a mathematical point of view, it is convenient to regard all the 
small-disturbance theories as being asymptotic forms of the full theory 
_ - - - — — — 

If the slope is high in some small region of the body, as at a 
slightly blunt leading edge, it may be presumed that the small-disturbance 
theory remains valid except locally. In this case x might be taken 
to be the thickness ratio. 
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for vanishingly thin bodies. Thus, the criterion for hypersonic flow 


may be expressed more formally as : ^ 


iT;. = 0(1) as 


We introduce a Cartesian coordi- 
nate system with the positive x 
axis alined with the free -stream 
direction (sketch (a)). Let the 
surface of the body be described by 
B(x,y,z) = 0, and the complete system 
of shock waves by S(x,y, z) = 0, 
where the function S is not, of 
course, known at the outset. 


{ - ( 2 ) 



Sketch (a).- Notation for hyper- 
sonic flow past thin body. 


Full Problem 


Consider the problem of determing the three velocity components 
V, w, pressure p, and density p throughout the flow field in the 
vicinity of the body. The mathematical system required is the differen- 
tial equations of motion (which govern the flow except at discontinuities), 
the Rankine-Hugonlot relations across shock discontinuities, and boundary 
conditions at the surface of the body and far from the body. 

Equations of motion .- The differential equations of motion, which 
express the principles of conservation of mass, momentum, and energy, 
are 

(continuity) (pu)^ + (pv )y + (pw)^ = 0 (3a) 

(x momentum) uu^^ + vuy + WU 2 + Px/P = 0 (3b) 


^Here the order symbols are used in the conventional sense: 
f(r) = 0 ( 1 ) as T— >0 means that f(r) remains bounded as r—> 0; 
and f(T) = 0(g[r]) means f(T)/g(r) is 0(l); similarly, f(r) = o(l) 
means f(T) vanishes as t — > 0. 
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(y momentum) 

UVx + Wy + WVg + Py/p = 0 

(3c) 

(z momentum) 

uwx + vWy + ww^ + p^/p = 0 

(3d) 

(energy) 

u(p/p’')^ + v(p/p^)y + w(p/p^)^ = 0 

(3e) 


(See, e.g., ref. 9 , ch. 1.) Here subscripts Indicate differentiation, 
and 7 is the adiabatic exponent of the gas. The last equation actually 
expresses the fact that the entropy is constant along streamlines, which 
for steady flow is equivalent to the conservation of energy (ref. 9 } 

pp. 15-16). 

Boundary conditions . At the body the normal component of velocity 
must vanish. The unit normal vector at the surface is proportional to 
grad B, so that the condition becomes 

q • grad B = 0 


where q is the velocity vector, or 


(tangency) uBx + vBy + wBg =0atB=0 ( 4 ) 


The other boundary condition, which implies that the body is flying 
into still air, may be taken in various equivalent forms. For present 
purposes it is convenient to require that all disturbances vanish far 
ahead of the body: 


u u. 


(upstream) v and w ->0 

w 

p -> p_ 


\ as 


( 5 ) 


Shock-wave relations .- At a shock wave, conservation of tangential 
momentum leads to the requirement that the velocity component tangent 
to the shock surface be continuous. The tangential velocity component 

(sketch (b)) is given by q^ = (n x q) x n, where the \init normal 
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vector n is proportional to grad S. 
It is convenient to use brackets to 
denote the jump in a flow quantity 
across a shock wave so that, for 
example, [u] is the increase in u 
through the shock. With this nota- 
tion, the condition of conservation 
of tangential momentum becomes 

(tangential momentum) 

M = M = M at S = 0 ( 6 a) 

Sx Sy Sg 


This imposes two independent scalar conditions at the shock wave, as 
physical considerations clearly indicate that it should. For plane 
flow, say, in the x-y plane, the last term becomes indeterminate and 
should, of course, be dropped. 

The remaining shock-wave relations express the conservation of nor- 
mal momentum, mass, and energy across the shock. The magnitude of the 
velocity component normal to the shock is (sketch (b)). 


y 



Sketch (b).- Components of 
velocity at shock wave. 


^ ^ uS +vS^+wS_ 

In = <1'^ = ^ 

ySx^+Sy^+Sz^ 


Consequently, the other three shock relations are found to be (ref. 9 ^ 
p. 300), using the jump notation. 


(mass) [p(uSx+vSy+wS2)] = 0 




( 6 b) 


(normal 

momentum) 


(energy) 


[p(uS^+vSy+wS2)^ + (S^^+Sy^+S2^)p] = ® 
|(uSx+vSy+wS2)^ + ^(S/+Sy2+S2^)£ 


Sat S =0 



( 6 c) 

( 6 d) 


So far, these relations are quite symmetrical, remaining unchanged if 
the brackets are taken to denote the change upstream rather than down- 
stream through the shock. A definite sense of flow direction is 
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provided only by the second law of thermodynamics, which requires that 
the entropy shall not decrease across each shock wave, so that 


(2nd law of thermo.) >0 at S = 0 (6e) 

L P^- 


For later use we record Bernoulli's law 


1 (u^+v^+w^) + ^ = const. (7) 

2 7-1 P 

whxch holds, with the same constant, throughout the flow field. 


First-Order Problem 


The full problem is now to be simplified by discarding all but lead- 
ing terms in the body thickness t. This will give a first-order hyper- 
sonic small-dlst\irbance theory, which can be expected to provide a close 
approximation for thin shapes. 

The reduction is conveniently carried out by introducing new inde- 
pendent variables which are of order unity throughout the flow field. 

The form of this transformation is suggested by simple examples and 
limiting cases. For example, the approximate solution for a thin plane 
wedge at very high Mach number has been given several times. Pertinent 
results are that the lateral extent of the flow field is some moderate 
multiple of the wedge thickness, that the density never exceeds 
(7+l)/(7-l) times its free-stream value, and that the relative streamwise 
velocity disturbance and the pressvire coefficient are of the order of 
the square of the svirface slope. Again, the Newtonian impact theory, 
modified to include effects of centrifugal force, which yields the limit- 
ing solution for M = « and 7=1, shows pressure coefficients propor- 
tional to the sqxiare of the thickness. The tangency condition suggests 
that, in general, all cross-wind velocities vary directly with the thick- 
ness. Such considerations suggest introducting new (barred) independent 
and dependent variables, and redefining the functions describing the 
body and shock-wave surfaces, as follows: 
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u = 

V = 

w = 

p = 
p = 
B = 
S = 


The new dependent variables 
regarded as dimensionless if the body is of linit length. 

We tentatively assume that all the new dependent variables (u, v, etc . ) 
and the new fimctions B and S are 0(l) as x -> 0 for fixed Mx, and 
that the reciprocals of the new independent variables (l/x, etc.) are 
likewise 0(l). The correctness of this assumption is suggested by 
examples such as those discussed above; its justification will come 
from the consistency of the resulting theory. 

It is important to realize that the notation v = 0(l) includes 
the possibility that in the limit v becomes arbitrarily small as well 
as the possibility that it approaches a constant nonzero value; only the 
possibility of its growing arbitrarily large is ruled out. For example, 
the reduced velocity components u, v, and w will be identically zero 
in the region ahead of the body. On the other hand, it is definitely 
implied here that in at least some portion of the flow field the reduced 
quantities will not vanish in the limit as x -^O.^ To be sure, they 
jnay not all be of order of magnitude unity in the intuitive physical 
sense; for example, for flow past a thin flat wing u and v will be 
moderate multiples of unity, but w will be numerically much smaller. 

Reduced problem .- This transformation of variables is now intro- 
duced into the full problem of equations (3) to ( 6 ). If we discard 

terms which contain x^ explicitly, such as t^u^ in the continuity 
equation, the differential equations become 

(x momentum) u- + + wu- + p-/p = 0 (9) 


u«[l+T^ u(x,y,z)] 
u„x v(x,y,z) 
u„x w(x,y,z) 
P„(7M^x^) p(x,y,z) 

B(x,y,z) 

S(x,y,z) 




(8b) 


J 


(80) 


are dimensionless, and the others may be 


“^That is, f(r) = 0(x) implies that f(x) is not identically o(x). 



10 


NACA TN 3173 


(continuity) 

P- + (pv)- + (pw)- = 0 

(10a) 

(y momentian) 

V- + + w- + p- / p = 0 

X y z ^ 

(10b) 

(z momentum) 

w- + w- + w- + p-/d = 0 
X y z -^z^^ 

(10c) 

(energy) 

(p/p\ + v(p/pV -f w(p/p^) =: 0 

(lOd) 


the boundary conditions become 


(tangency) 

(upstream) 


B- + vB- + wB- = 0 at B = 0 
X y 2 


( 11 ) 


u, V, w 0 
p 

P 1 


) as X — > - » 


( 12 ) 


and the jump conditions at the shock waves become 


(tangential 

mcanentum) [u ] _ [v ] _ [w] 

§- §- ' §- 

X y 2 




(mass) [p (S- + + ^-) 1=0 

X y z 


(normal _ » _ _ _ 

momentum) [p(S- + vS:z + wS-)^ + (S-^ + S-^)P] = 0 

X V z y z 


(energy) 




at S=0 (13) 


i(S- + vS- + wS-)^ + (Sy^ + Sz^) £ 
2 ' X y z' 7-1 ^ z / p 


= 0 


(2nd law r p 
of thermo.) ® 


The parameters- M and t of the full problem enter this reduced 
problem only in the conibination Mt, which appears only in the upstream 
condition on p (eq. (12)). In the hypersonic range, where 1 /Mts=0(1), 
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the reduced problem possesses complete internal consistency. This is 
most readily vinderstood by considering first the case where it does not. 

In the range of linearized theory, where® Mr « 1, equation (12) 
presents a contradiction. It was assumed at the outset that p is 0(l), 
but its value upstream, (l/M^®), is then not 0(l). This inconsistency 
is an automatic warning that the reduction breaks down in the range of 
linearized theory. Indeed, it will be seen later that the assumed orders 
of reduced quantities are then actually incorrect, and that the reduced 
equations fail to describe linearized supersonic flow. 

In the hypersonic range no such inconsistencies arise. Hence, the 
tentative assumptions regarding orders of reduced quantities are justified 
a posteriori , as are the simplifications effected by discarding terms that 

involve explicitly. 

Order of error .- Because terms of order have been omitted in 

the reduction, first-order quantities will differ from their exact values 
by O(t^). For the special case of plane flow, an analysis similar to 
the preceding has been given by Goldsworthy (ref. 8), who considers only 
a single bow wave;® he furthermore confines attention to the range Mt ~l 
(which, to be sure, may be the range of most practical importance) and 
finds the error to be 0(1/M^). This is equivalent to the present result 
in that range. However, the present result is more general, holding for 
arbitrarily large values of the similarity parameter (assuming, of course, 
that the assumption of a continuimi flow remains valid). For example, at 
Infinite free-stream Mach number, the error in first-order theory is 
correctly 0 (t^ ) . 

It is interesting to note that the error in the various first-order 
small -disturbance theories decreases progressively from 0 (t^/®) in tran- 
sonic flow to 0(t) in linearized supersonic flow to O(t^) in hypersonic 
flow. Therefore, under the plausible assumption (confirmed by later 
examples) that these mathematical order estimates give a reasonable indi- 
cation of the actual physical magnitude of error, the practical need for 
a second-order solution is seen to be greatest for transonic flow and 
least at hypersonic speeds. 

Unsteady analogy .- A significant feature of the reduced problem is 
that the problem for the streamwise velocity u has been uncoupled from 
that for the other variables. Equations (lO) to (l3) constitute a com-_ 
plete problem for v, w, p, and p, which can be solved independent of u. 

®More precisely, where Mt = o(l) as t^O. 

® It seems important to consider multiple shocks, because it is not 
at all obvious that subsequent- shocks have the same status as the bow 
shock; in particular, at M = « the upstream Mach number is Infinite for 
the bow shock, but has moderate supersonic values for the subsequent 
shocks . 
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Thereafter, u (if required) can be determined from Bernoulli's law: 


u + ^ f = const. 

2 7-1 p 


(7-1)M^t^ 


(IM 


Consequently, equation (9) and the first terms of equations (l2) and 
(l3a) are superfluous and can henceforth be disregarded. 

As pointed out by Hayes (ref. 6), this first-order, problem is com- 
pletely equivalent to a full problem for unsteady flow in one less space 
dimension. The reduced problem of equations (lO) to (13) is precisely 
the full problem of unsteady motion in the y-z plane due to a moving 
solid boundary described by B(x,y,z) = 0, where x is interpreted as 
the time, and all other barred variables as the actual physical quan- 
tities."^ The outline of the moving boundary is given by the trace of 
the original thin shape in a cross-stream plane which moves downstream 
with the free -stream velocity (sketch (c)). For example, the problem 
of steady hypersonic flow past a slender pointed body of revolution is 

equivalent to the problem of unsteady 
planar motion due to a circular 
cylinder whose radius varies with 
time', growing from zero at time 
X = 0. Hayes has given a physical 
explanation of this analogy. 

It may be noted that the analogy 
is similar to that arising in the 
slender-body theory of linearized 
compressible flow, as exemplified by 
the work of Jones (ref. 10 ) and Ward 
(ref. ll) . There, however, the time- 
dependent analogue is incompressible, 
whereas here it is definitely com- 
pressible. 


Similarity rule .- The parameters M and t appear in the reduced 
problem only in the combination Mt, which is the hypersonic similarity 
parameter of Tsien. This means that for bodies derivable from one 
another by uniform contraction or expansion of all dimensions normal to 
the stream, the flow fields are related if the corresponding Mach 



Sketch (c).- Plane of . 
unsteady analogy. 


^The undisturbed fluid has density pressure p^=i/ 7 M^t^. 
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numbers are such that the similarity parameter Mr is the same in each 
case. The nature of this relationship is simply that the reduced flow 
quantities, as functions of the reduced coordinates, are the same in each 
flow. All distinctive flow surfaces, such as the body itself and the 
shock-wave system, have identical descriptions in terms of the reduced 
coordinates. This is the hypersonic similarity rule of Tsien, which was 
extended to rotational and three-dimensional flow by Hayes. Its implica- 
tions with regard to pressure and force coefficients will be summarized 
later, after the hypersonic rule has been combined with that for ordinary 
supersonic flow. 

Extension to unsteady flow .- For simplicity, the preceding discussion 
has been restricted to steady flow. It can readily be extended to unsteady 
motions involving small time-dependent oscillations of a thin body exposed 
to a steady uniform stream (or, from another point of view, flying through 
still air executing sli^t time -dependent variations from a mean steady 
rectilinear fli^t). The full problem is obtained by replacing the sub- 
stantial derivative u^/dx + vb/dy + vb/bz wherever it appears in equa- 
tions (3) to ( 6 ) by its unsteady counterpart b/bt + \ib/bx + vb/by + wb/bz. 
Correspondingly, if a reduced time t is introduced according to 

t = u„t (15) 


then the unsteady smal 1 -disturbance problem is obtained from the steady 
problem of equations (9) to (l 3 ) by adding b/bi throughout to the 
operator b/bx + vb/by + wb/bz. For example, the reduced continuity 
equation (eq. (lOa)) becomes, for unsteady motion 

+ % + (P^)y + (P^)g = 0 (16) 


The problem for u remains uncoupled from the problem for _ the other 
variables.® Fxirthermore , in the remaining problem the x and t deriv- 
atives appear only in the combination d/dt + b/bx. This means that the 
number of independent variables can be reduced by one by adopting the 
viewpoint of an observer who is moving with the free stream (or, to put 
it another way, who is fixed in still air with the body moving past). 
Thus, introducing "still-air" coordinates 


X ' = X - u*t X ' = X - u„t 

t' = t t* = t 


(17) 


®However, u can no longer be found immediately in terms of the 
solution of the remai n ing problem, because in unsteady motion there is 
no useful counterpart of the Bernoiilli equation (eq. ( 1 ^ 4 -)); fortvinately, 
u is seldom actually required. 


/ 
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reduces the unsteady small -disturbance problem to exactly the form of 
the steady problem (eqs. ( 9 ) to (13) with t replaced by t'. This 
means that Hayes ' analogy remains valid for unsteady motion if account 
is taken of the actual variation of the contour with time as well as the 
apparent variation due to relative motion depicted in sketch (c). This 
result was first given by Hamaker and Wong (ref. 7)» Recently, Lighthill 
has analyzed oscillating airfoils at hypersonic speeds from this point of 
view (ref. 12). 


If the body oscillates so rapidly or with such large amplitude that 

B- is much greater than B-, the error remains of O(t^) only if t is 
t ^ 

taken to be the maximum Instantaneous slope of particle paths at the sur- 
face of the body. 


Unified supersonic-hypersonic small-disturbance theory .- It would 
be advantageous if the hypersonic small-disturbance theory included 
linearized supersonic theory as a special case. Then the awkward ques- 
tion of what is the lower limit of hypersonic flow would not arise. 

A single theory, and corresponding similarity rule, would hold for all 
supersonic speeds above the transonic range. 

In the case of transonic small-disturbance theory, a unification of 
this sort is known to arise quite naturally (ref. 13 , p. 9)« In its 
original form (with not replaced by 2(M-1)) transonic theory 

embraces linearized theory as a special case, so that it furnishes 
effectively a unified subsonic-transonic-supersonic theory giving a 
first approximation at all speeds below the hypersonic range. 

In the case of hypersonic flow, a connection with the adjoining 
supersonic range can likewise be effected (ref. 1^), but the reason 
therefor is much less strai^t forward. An immediate obstacle is the 
fact that the approximations leading to the hypersonic theory and those 
leading to linearized theory are mutually exclusive. The difficulty 
arises in the continuity equation (eq. (3a)) which in the hypersonic 
case was shown to reduce to 

u p + ( pv) + ( pw) =0 ( 18 a ) 

but in linearized theory reduces instead to 

^coRx + + Pj’^y + ^z) = ° 

The term p^Uj^ must be retained in linearized theory, whereas it must 
be neglected in hypersonic theory in order to achieve similitude. It 
would therefore appear impossible to give a small-distiorbance theory 
which is general enough to describe both hypersonic and ordinary super- 
sonic flows, and yet simplified enough to retain the corresponding 
similitudes. It must be regarded as a coincidence that this can, 
nevertheless, be accomplished. 
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It is shown in Appendix B that the hypersonic theory covers the 
ordinary supersonic range if it is reinterpreted in accordance with the 
similarity rule for linearized theory (ref. 15). It is found that solu- 
tions of hypersonic small-disturbance theory remain valid at small 
values of the parameter Mr (which is the domain of linearized theory), 
provided that Mr is replaced by and the results are reinterpreted 

in terms of physical variables according to 


u = u^[l + u(x,y,S; pr)] 
V = u«T v(x,y,z; pr) 


w = u«T ^?(x,y,z; pr) 


P = P» p(x,y,Z; pT) - — ^ I = p^^7MV 

P = + 3t) - 1 ]} = -i) 



rather than according to equation (8). The pressure coefficient is 
given by 





(19b) 


Since the error is O(t^) in the hypersonic theory and 0(x /p ) in the 
linearized theory, the error in this unified theory is O(t^) or 0 (t/P), 
whichever is the greater.® 

The flow quantities of chief aerodynamic interest are the pressure 
coefficient and the various force and moment coefficients derived from 
it by integration. If the hypersonic small -disturbance problem has 
been solved to find Cp/r^ as a function of the reduced coordinates 
and the parameter Mr, the result is rendered valid also in the ordinary 
supersonic range simply by replacing Mr by pr . 


Unified similarity rule .- The unified supersonic-hypersonic similar- 
ity rule may be summarized as follows: 

For steady flow past thin bodies derivable from one another by uni- 
form contraction or expansion of all dimensions normal to the stream, 

®For slender shapes such as th: 
is only 0 (t^/P^) in the linearized 
0 (t^/P^) in the combined theory. 


n bodies of revolution, the error 
theory and, hence, O(t^) or 
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the flow fields are related if the corresponding Mach mmbers are such 
that the similarity parameter Pt is the same in each case (t being 
any measure of thickness). The relationship is such that the flow 
fields are identical when expressed in terms of the reduced flow quan- 
tities 


u-u^ V w p-p^ 

i ) ) p p 

T “^U^ TU„ TU„ M‘=^T‘='p^ 

as functions of the reduced coordinates 


When the body oscillates slightly, the same is true with the addition of 
u^t to the reduced coordinates, provided that the time history of oscil- 
lation in terms of reduced coordinates is the same for each body. 

From the rule for pressures follows the equality of the reduced force 
and moment coefficients 



2+k* _i+k^ * 

T T T 


etc . 


Here k = 1 if some plan-form area is taken for reference and k = 0 
if some cross-sectional area is used. The connection between the simi- 
larity rules for forces and moments is contained in the statement that 
the center of pressure is constant in terms of the reduced coordinates. 


An vmlimited number of equivalent forms of the reduced variables 
can be produced, for example, by multiplying by any power of the simi- 
larity parameter. However, the forms given here are the most useful ones 
in the hypersonic range, because they involve functions of order unity. 

In the supersonic range no forms have this advantage, except in the 
special case of plane flow, where it would be convenient to use 



It should be noted that in the hypersonic and combined similarity 
rules the adiabatic exponent y must remain fixed, whereas its magni- 
tude is arbitrary in the supersonic case (since it does not appear in 
the linearized problem). 
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TYPICAL APPLICATIONS OF HYPERSONIC 
SMALL-DISTURBANCE THEORY 


Several problems will now be solved according to hypersonic small- 
disturbance theory. These examples will illustrate possible methods of 
solution, and demonstrate the degree of simplification resulting from 
the assumption of small disturbances. Comparison with the corresponding 
solutions of the full equations (when available) will Indicate the 
accuracy to be anticipated when the theory is applied to more elaborate 
problems . 


Plane and Axially Symmetric Flows 


The examples to be considered are either plane or axially symmetric 
flows. Accordingly, it is convenient to introduce coordinates x,r 
where in the case of plane flow r is the Cartesian coordinate y. 
Henceforth, v will denote the velocity component in the r direction, 
which is the radial velocity for axially symmetric flows. 

Equations of motion .- In these coordinates, the hypersonic small- 
disturbance equations of motion (eqs. (lO),) become 


+ (pv)j + a ^ = 0 (20a) 

V- + w- + — = 0 (20b) 

X r ^ 

(p/p^)- + v(p/p^)j = 0 (20c) 

The distinction between two and three dimensions arises only in the 
continuity equation, where a = 0 for plane flow and a = 1 for 
axially symmetric flow. 

Boundary and shock conditions. - 
Let the surface of the body be 
described by r = Tb(x), as indi- 
cated in sketch (d). Then the 
tangency condition of equation (ll) 
becomes 

V = b*(x) at r = b(x) (21) 

The upstream conditions are given 
by equation (12). 
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The examples to be considered will involve only a single bow shock 
wave, which may be described by r = Ts(x). The shock-wave conditions 
are given by equations (13), but in this case it is easier to take advan- 
tage of the convenient relations which are listed in reference 16. 
Equations (136), (128), and (129) of reference I6, when reduced to hyper- 
sonic small-disturbance form, give 


V 


2 K ^-1 

7+1 


s'(x) 


p = 


- (y-1) 

r(7+i)K^ 


^(x) 




at f = s(x) 


0 = (7+1) ^ 

2 + (7-1)k® 

J 


( 22 ) 


where K = Mts*(x) is the hypersonic similarity parameter based on the 
local shock-wave slope. 


Stream function .- The continuity equation (eq. (20a)) may be 
accounted for by introducing a stream function, setting 



a 

r pv = -+£ 


(23) 


so that 



(24) 


.Then equation (20c) states that p/p' is a function o nly of i|f, as is 
clear from the fact that entropy is constant along streamlines between 
shock waves. Substituting into the momentum equation (20b) gives 


7+1 

1 r- (7wt” + 


- 2Vf%+V'^ff =< 


for plane 
flow 


( 25 a) 




for 

axially ( 25 b) 

symmetric 

flow 
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where 


w(t) = p/p’^ 

The pressure is given by 


(25c) 


( 26 ) 


Plane Wedge 


As a simple introductory example, consider hypersonic flow past a 
thin wedge of semi vertex angle 6. Here, and in the examples to follow, 
the solution is most readily carried out by assuming a given shock wave 
and Mach nxanber and calculating the corresponding body shape. It is 
therefore convenient to identify the thickness parameter T with the 
shock-wave angle rather them the wedge angle. This is quite permissible 
at hypersonic speeds, where they eire of the same order; it is also per- 
missible in the ordinary supersonic range, where they are not, provided 
that 8 rather than t is used in the error estimates. Let b = 8/t 
be the ratio of wedge to shock angles. 

The flow field is conical (velocities constant along rays), so that 
the stream fxmction has the form 


i|r(x,r) = X f(e) 
where 6 is a conical variable, defined by 



(27) 

( 28 ) 


which varies from b at the wedge to unity at the shock wave . Thus , 
equation (25a) becomes 


f" 




= 0 


(29a) 


where, from equations (22) and (25c) 


“o = 


27k/ - (7-1) 


7(7+1) Ko^ L (7+1) Ko^ 


2 + (7-1) ^7 


(29b) 


with Kq = Mt. The shock-wave conditions of equation (22) give (since 
s* = 1) 

f(l) = 1 

(7+1) K ^ } (30) 

f»(i) = — L.0 — f > 

2 + (7-1) Ko^ 
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The solution of this problem, which corresponds to f" = 0, is 


f(e) 


(y+1) k/ e - 2(Kq 2 -1) 
2 + (7-1) Ko^ 


(31) 


Requiring f to vanish at the surface gives the ratio of wedge to shock 
angles 


8 _ ^^*^0 ~ 

T (7+1) Ko^ 


(32) 


The auxiliary hypersonic parameter Kq = Mr can now be eliminated in 
favor of M6, with the result that 



(33) 


Then from equations (26) and (27) the pressure coefficient at the wedge 
(or anywhere between the wedge and shock wave) is found to be 



(3i^) 



Sketch (e).- Wedge pressure 
according to hypersonic 
small-disturbance theory. 


These results were first given by 
Linnell (ref.lT)* Numerical values 
of pressure coefficient are com- 
pared in sketch (e) with the exact 
results for wedges of various thick- 
nesses. Here 6 has been taken to 
be the tangent of the wedge angle, 
though within the scope of the 
small-disturbance theory it might 
equally well have been Identified 
with the sine of the angle, the 
angle itself, etc. 
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Unified supersonic-hypersonic result .- Replacing the parameter M6 
by 36 renders the solution valid in the ordinary supersonic range as 


well: 



(35) 


This result is again compared with 
the full solutions in sketch (f). 
The advantage of the unified resvilt 
is obvious. In all subsequent 
examples the results will be pre- 
sented only in this form. 

It is interesting to note that 
this formula has been proposed by 
Ivey and Cline (ref. I8) for pre- 
dicting the surface pressure on any 
supersonic airfoil. They obtained 
it by seeking an interpolation 
formula connecting Ackeret's linear- 
ized theory with the hypersonic 
result of Linnell (eq. (3^)) for a 
tangent wedge. 



Sketch (f).- Wedge press\ire accord- 
ing to unified supersonic - 
hypersonic small-disturbance 
theory. 


Initial Gradients for Plane Ogive 


Using the full equations, Crocco first determined the initial gra- 
dients of flow quantities at the tip of a plane ogive by perturbing the 
solution for a wedge (ref. 19). His analysis has been repeated and 
elaborated upon by Schafer (ref. 20) and others. This problem provides 
a good test of the small -disturbance theory in a case involving shock- 
wave curvature. Let b be the ini- 


tial ratio of body slope to shock- 
wave slope, and I the correspond- 
ing ratio of radii of cixrvature. 
Then in physical coordinates the 
body (sketch (g)) may be described 

by 

r^ = r(bx + |cx^ + ...) (36) 

and the shock wave by 



rs = t(x + |zcx^ + ...) (37) 


Sketch (g).- Plane or axially 
symmetric ogive. 
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Here c will be negative for the convex shapes usually encountered in 
practice. 


Shock-wave conditions .- Conditions just behind the shock wave are 
found from equations (22) to be given by 



pv 


2(Kq^- 1) f ^ (7-1 )Kq-^+(7+5)Ko^- 2 

2+(7-1)ko^ ^ [2+(7-1)Ko^] 


lex + 



( 38 a) 


(7+1)Kq2 

2+(7-1)Kq^ 


1 + ^ lex + . . . 

L 2+(7-1)Kq2 


( 38 b) 



27I«o^-(7-1) 

7(7+1)Kq^ 


2+(7-1)Kq2 

(7+1)Kq^ 



^ 7 ( 7 - 1 ) (Kq^- 1 )^ 

lex + 

[27Ko^-(7-1)] [2+(7-1)ko^ ] 


} 


(38c) 


where Kq = Mr is the auxiliary hypersonic similarity parameter based 
upon the initial shock-wave angle. 

Equations of motion .- The flow behind the shock wave will consist 
of a uniform field upon which is superimposed a perturbation field due 
to body • curvature . Hence, along each ray from the vertex the flow quan- 
tities will have constant values associated with the initial slope of 
the body plus linear variations proportional to the initial curvat\ire 
(together with higher variations which need not be considered in evalu- 
ating initial gradients). Thus, the stream function of equation (25a) 
can be written in the form 

i{f(x,r) = X f(e) - Icx^ g(e) + ... ( 39 ) 


Along the shock wave 


t = f = X + 


• • « 
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so that the entropy function w can be written as ilr(w) by replacing 
X by in equation ( 38 c). ■ In this fonn (since U is constant along 
streamlines) the expression applies throughout the flow field downstream 
of the shock wave, so that 

^ c 

= Wq[ 1 + Icuiilf + ...] _ (ij-Oa) 


where 


Wr 


27Kq2-(7-1) 

7(7+1)Kq^ 


~ 2+(7-1)Kq^ ~ 
. (7+1)kq^ _ 


7 


^ 7 ( 7 - 1 ) (Kq^-1)^ 

[27Ko^-(7-1)1 [2+(7-1)Ko^] 






) 


(il-Ob) 


Substituting these expansions for ^ and 0 ) into equation ( 25 a) 
and equating like powers of x yields two ordinary differential equa- 
tions. The first is equation (29a), corresponding to the basic flow 
past a wedge. The second, when simplified with the aid of equation (31), 
becomes 

fV-2ff'g'+2f’2g = (7g"-Uif'2) (4l) 


Solution for g .- By expressing conditions just behind the shock 
wave in terms of conditions at 0=1 through Taylor expansion, the 
boundary conditions on g are found to be ^ 


g(l) 


2+(7-1)i«o2 


g'(l) = - 


lf(7+l)K/ 

[2+(7-1)'«o^ f 




J 


(42) 


The velocities associated with g vary linearly with distance from 
the apex. Hence (as is readily verified from the differential equation), 
g is a quadratic function of 0 , and can therefore be written as 


g(0) = g(l) - (l-0)g’(l) + |(1-0)V(1) (i^3) 


Here only g”(l) is unknown, and it is immediately found from the ‘dif- 
ferential equation to be 


g”(l) = 2(7+1 )ko^ 


3 ( 7-1 ) + ( 3-7 ) Kq^+2 ( 7+2 ) 

(Ko^-l) [2+(7-1 )Ko^]® 
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Curvatiire ratio and siirface pressitre .- The body curvature and pres- 
sure gradient can be expressed in terms of the values of g(d) and its 
derivative at 0 = b. From equations (42) to (44) it is found that 


■ ' g(b) 


2(27-1)Ko^+(7+^)Kq^-(7-1) 
(7+1)(Ko^-1) [2+(7-1)Ko^] 


s' (b) 


2k/(3Kq^-H) 
(Kq^.i) [2+(7-1)k/] 


> 


J 




The surface of the body is determined by the vanishing of the 
stream function. Thus it is found that the ratio of shock-wave curvature 
to body curvature is given by 


I = 

2g(b) 


(7-H)^k/(k/-1) 

2 [2(27-1)Ko‘^+(7+5)Ko^-(7-1) ] 


(46) 


The initial pressure gradient on the surface of the body is given 


by 


* = -7(»),f'(b)^ g'(b) T^c 

g(b) 


■^P 

6x 


(^7) 


so that in terms of the initial slope Rq' = br and curvature Rq" = ct 
of the body 


6c 


P 

6x 


(7+1)i«o^(3Ko^+1) [27Ko^-(7-1)] 
(Kq^.i) [2(27-1)Ko^+(7+5)Ko^-(7-1)1 


- Ro'^o" 


(48) 


Although these results have been expressed in terms of the auxiliary 
parameter = Mr, they can be given explicitly in terms of Mach number 
and apex angle (in the combination MRq' = Mbr) with the aid of equa- 
tion (33). In this respect the small-disturbance solution is superior 
to the full solution, which yields no such explicit results. Replacing 
MRq’ by PRq' renders these results applicable at all supersonic speeds. 
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The small-distiorbance result for the curvature ratio is compared 
in sketch (h) with the full solution for various vertex angles, taken 
from the convenient tabulation of reference 21 . Sketch (i) gives the 
corresponding comparison for the initial pressure gradient. 



Sketch (h).- Initial ratio of 
shock to body curvature 
for plane ogive. 



Sketch (i).- Initial pressirre 
gradient on plane ogive • 


Circular Cone 


Consider flow past a slender 
circular cone of semivertex angle 6 , 
(sketch (j)). Again it is conven- 
ient to regard the Mach nmber and 
shock-wave angle as given, and to 
solve for the corresponding cone 
eingle. Thus let the shock-wave 
angle be x and the cone angle 
6 = hr, where again b is a con- 
stant less than unity that is to be 
determined. 



Sketch (j).- Notation for cone. 


Equation of motion .- The flow field is conical so that all flow 
quantities are constant along rays. This means that the stream function 
of equation (25b) has the form 


t(x,f) = x 2 f(e). 





It follows that the equation of motion 
ordinary differential equation 

i^fSfti _ 2ff'^ = ro)o — 

0 


(eq. 

(7+x) 

(7-1) 


(25b ) ) becomes the nonlinear 

0" - V) 


(50) 
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Here % is the constant value of p/p^ behind the shock wave, given 
by equation (4ob), and Kq = Mr is again the auxiliary hypersonic simi- 
larity parameter based upon shock-wave angle. 


Boundary- conditions . - The shock-wave conditions of equation (22), 
together with equation (U9) , combine to give 


f(l) 



f’(l) 


(y+l)Ho^ 

2+(7-1)k/ 


(51) 


The condition of tangent flow at the cone requires that the stream 
function vanish at the surface: 


f(b) = 0 (52) 

Numerical integration ,- The nonlinear equation for equation (50)^ 
can be readily integrated numerically. Choosing a value of the auxiliary 
similarity parameter we calculate o)q from equation (40b) and the 

initial values of f and f* from equations (51)^ and then integrate 
step by step inward from the shock wave until f vanishes^ which deter- 
mines the cone surface. With the ratio of cone angle to shock-wave 
angle 6/t = b thus determined^ the results can be re-expressed in 
terms of the similarity parameter based upon cone angle. Eight or ten 
intervals between shock wave and body yield ample numerical accuracy^ 
provided that in each step the predicted values of f and its deriva- 
tives are corrected by averaging and iterating before proceeding to the 
next step. (For values of Kq near unity ^ the first few intervals 
near the shock wave must be taken smaller than the others.) 

The pressirre coefficient is obtained in terms of the first deriva- 
tive of the stream function according to 


^ = 2 

r“o 

f—') 

1 

H 

1 

6^ 

_b2* 

\e J 

J 


( 53 ) 


where the similarity parameter M6 of the hypersonic problem is to be 
reinterpreted as p& so that the result is applicable throughout the 
entire supersonic range. 

Computations have been carried out for y = 1.405, in order to com- 
pare with the full solutions tabulated by Kopal (ref. 22). The chosen 
values of the parameters are listed in the following table, together 
with certain of the results : 
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Mt 

6 / T 

p6 

f(b)/b 

CM 

0 

1.04 

1.19 

1.58 

2.87 

4.47 

00 

0.3620 

.5545 

.7281 

. 86 o 4 

.8921 

.9i4o 

0.3765 

.6599 

1.150 

2.469 

3.988 

00 

1.217 

1.522 

2.207 

3.925 

4.973 

6.149 

3.183 

2.646 

2.333 

2.154 

2.116 

2.091 


The small-disturbance result for surface pressure is compared with 
the full results (from ref. 22) in sketch (k). The differences between 


the full solution and the small - 
disturbance limit are closely pro- 
portional to the square of the 
thickness^ in accordance with the 
estimate of the error as 0(&^) 
or 0(6^/p^). It is noteworthy that 
the fractional differences are in 
fact very nearly equal to 6^. 

The same is true of surface pres- 
sures in the previous examples. 

This suggests that the mathematical 
order estimate may be relied upon 
to give a quantitative prediction 
of the error in the small- 
disturbance approximation. 

This approximate solution for 
cones was previously sought by 
Shen (ref. 23 ), whose result is 
also shown in sketch (k). It 
appears that his solution, which 
involves more Involved computa- 
tions, must contain errors. 



Sketch (k).- Surface pressure 
on cone . 


Initial Grad.ients on Ogive of Revolution 


Consider the axially symmetric counterpart of Crocco’s problem: 
the determination of the initial gradients of flow quantities at the 
tip of an ogival body of revolution. This problem has recently been 
considered for the full equations by Cabannes (ref. 24) and by Shen 
and Lin (ref. 25). It will be seen that the small -disturbance solution 
serves as a useful guide to the full solution. In particular, it clar- 
ifies the behavior of the solution near the surface, which was glossed 
over by Cabannes and in Shen and Lin's work was misconstrued to 
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indicate a singularity which implied wrongly that the initial pressure 
gradient at the tip of an ogive is infinite. 

As in the plane problem, let the body be described in physical 
coordinates by 

r^j = r(bx + i cx^ + ...) (5^) 

and the shock wave by 

rg = t(x + I Zcx^ ...) (55) 

(sketch (g)). Conditions just behind the shock wave are given by equa- 
tion (38), with -ijr- and replaced by -^-/r and t-/f. 

XI* X i 

Equations of motion .- Again along each ray the flow quantities 
have constant values corresponding to the initial slope plus linear 
variations proportional to the initial curvature (together with higher 
variations which are considered later)* Hence ^ the stream function may 
be written in the form 

\lf(x,f) = x^f(e) - cx°g(e) + ... (56) 

Along the shock wave 

= if ^ = ix^ + . . . 

2 2 

so that an expression for the entropy function u throughout the flow 
field downstream of the shock wave is obtained by replacing x with 
^ 2^ in equation ( 38c ) . Hence , 


w(\|/) = Uq [1 + ZcWi + ...] (57) 

J 

where Wq and cji are given by equation (40b). 

Substituting these expansions into equation (25b) and equating like 
powers of x yields for f the nonlinear ordinary differential equa- 
tion already treated in the problem of the cone (eq. (50))^ for g 

the linear equation 

Dg" = A + Bg + Cg' (58a) 


whose coefficients depend upon f according to 
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A = (OqU)! 


f ’ 


(7+1) 


e(y-^) 


- 

_/2f 





B = 12ff” 


C 


= 7Wo 


f,(r+i) 

0(7-1) 


7+2 

0 



8ff 


D = 7^0 


ft (7+1) 

e(r-i) 





(58b) 


As in the plane problem^ initial conditions on g are found by 
expressing conditions just behind the shock wave in terms of conditions 
at 0=1 through Taylor series expansions : 


g(i) = 




2+(7-1)k. 


g«(l) = - K( 


( 7-1 )Kq^+ 2 (27+3) 

[2+(7-1)k 2]2 


(59) 


Behavior of solution near body surface .- Just as conditions at the 
curved shock wave have been related to those at 0=1, so with the 
present coordinate system it is necessary to relate conditions at the 
surface of the body to those on the initially tangent cone 0 = b • 
However, the solution is nonanalytic near the surface, which means that 
Taylor series expansions do not exist. It is therefore necessary to 
examine the nature of the solution in the vicinity of the body. 

The fimction f associated with the basic conical flow is analytic 
near the surface (and vanishes at 0 = b), so that it has a Taylor series 
expansion : 

f(e) = (0-b) f(b) + ... O[(0-b)2] (60) 

It follows from equations (58b) and (60) that near the surface the 
coefficients of the differential equation for g behave like 

A ~ (0-b)"^/^ 

B ~ (0-b) 

C 1 f 

D ~ 1 


(61) 
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Therefore the point 0 = b is a singular point of equation (58a), but 
is an ordinary (nonsingular) point of the homogeneous equation obtained 
by deleting A (see ref. 26 , p. 73 ) • Therefore the general solution of 
the homogeneous equation is analytic and can be taken to have the form 



(62) 


where all higher coefficients an can be expressed in terms of the two 
arbitrary constants ao and a^ by means of the differential equation. 
Then the procedure for calculating a particular integral (ref* 26 ^ 
p* 122) shows that the nonhomogeneous equation has a particular 

integral of the form 


gp = ^ Cn(e-b)'' (63) 

o 

where the coefficients c^ can all be determined. Here the 3 / 2 -power 
branch point arises from the fact that the pencil of fluid striking the 
tip of the ogive is spread thin over the entire surface, and the linear 
entropy gradient at the tip due to a curved shock wave is thus intensi- 
fied to a square-root gradient normal to the surface elsewhere. The 
complete solution of equation (58) is the sum of gj^ and gp. 

In treating the full problem in reference 25, Shen and Lin claimed 
to have found a logarithmic singularity at 0 = b, which considerably 
complicated their analysis. Because of this singularity, their solu- 
tion was restricted to concave bodies (although they conjectured that 
it might be extended to convex bodies). The singularity also led to 
the conclusion that for an analytic body shape the shock wave is non- 
analytic, and vice versa. Furthermore, the singularity would .imply that 
the initial pressure gradient on an analytic body of revolution is 
infinite, although numerical solutions by the method of characteristics 
give no indication of this. 

The present solution shows no such singularity. It seems unlikely 
that a singularity could have disappeared as a result of making the 
small-disturbance approximation, since this would imply that the approx- 
imate model does not retain the essential features of the full problem. 
The alternative conclusion is that the singularity does not actually 
exist in the full problem. This has subsequently been confirmed by the 
authors of reference 25^° who find that the singularity is, in fact, 
only apparent in the sense of reference 26, page i|-06. 
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In a private communication; see also Addendum No. 1 to ref. 25. 



NACA TN 3173 


31 


It may be noted that Cabannes, in treating the full problem, has 
completely ignored the nonanalytic nature of the solution, and simply 
extrapolated his numerical solution to the surface of the body (ref. 24). 
His results may not be seriously in error, because the effect of the 
nonanalyticity is small. 

Numerical integration .- The differential equation for g has been 
integrated numerically for the six values of the similarity parameter 
chosen previously for the cone. The Integration was carried out step by 
step starting from the known values at the shock wave and using the same 
intervals as for the cone. This step-by-step solution was joined at the 
two points nearest the surface with the series expansion about 0 = b 
given by equations ( 62 ) and ( 63 ): 


g(e) 







g(b) [1 + 2 Xe^ - A(l+ 2 A)e'^ + . . .] + 


where 


g’(b) ( 0 -b) 




3-7 

~T 


Ae® + 



A = 



f*(b) 


(7-1)1- 1 


(7+1) J 


e-b 

b 


(64a) 

(64b) 


(Values from the step-by-step integration and the series expansion were 
also compared at the third point from the surface as a check.) Because 
they are based upon the previous solution for a cone, the computations 
were carried out \/ith 7 = 1.405. 


Curvature ratio and pressure gradient .- Because the nonanalyticity 
appears only in higher terms of the series, surface values of g and 
its first derivative (but no higher derivatives) can be expressed in 
terms of values at 0 = b. The surface of the body is determined by the 
vanishing of the stream function. Thus the ratio of shock- wave curvature 
to body curvature is found to be given by 


= f'(b) 

2 g(b) 


( 65 ) 


Proceeding as in the plane case, it is found that the initial pressure 
gradient is expressed in terms of the initial slope Rq' = bx and curv- 
ature Rq" = c t of the body by 
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bx 


f'(b)' 


M‘='Rr 


(7-3) 


[ -g»(b)l 
g(b) 


Ro’Ro" 


( 66 ) 


Numerical values of g(b) and g'(b) are listed in the following 
table ^ together with the resulting values for the curvature ratio I 
and surface pressure gradient: 


PRo' 

g(b) 

g' (b) 

l 

1 1 

Ro'Rq" 8 x 

0.3765 

5.170 

-8.502 

0.0426 

5.532 

.6599 

1.800 

- 3.640 

.2346 

4.662 

1.150 

1.573 

-4.073 

.5106 

4.514 

2.1^69 

2.101 

-6.820 

.8039 

4.684 

3.988 

2.487 

-8.638 

.8921 

4.802 


2.931 

-10.76 

.9586 

4.929 


The curvature ratio is plotted in sketch (z),and the initial pres- 
sure gradient in sketch (m) • The. curvature ratios calculated from the 



Sketch (i).- Initial ratio of 
shock to body curvature for 
ogive of revolution. 



Sketch (m).- Initial pressirre 
gradient on ogive of revo- 
lution. 


full equations by Shen and Lin are also shown in sketch (Z) for com- 
parison, because the error introduced by incorrect treatment' of the 
solution near the body is probably small. 
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Higher Terms in Series for Ogive of Revolution 


It might he supposed that for an ogive of revolution described by a 
power series, the perturbation scheme could be continued indefinitely to 
find successively higher terms in a power series expansion for surface 
pressure. However, because of the nonanalytic nature of the stream func- 
tion near the body, complications arise if one proceeds simply by adding 
f\irther terms to equation (56). 

It has been seen that g(0) involves an authentic 3/2-power branch 
point near the s\irface. However, it can be shown that the next term will 
involve a spurious l/2-power branch point there, the next an inverse 1/2- 
power singularity, and so on. As a consequence of this spurious reinforce- 
ment of the actual nonanalyticity, it is impossible to evaluate surface 
pressures. Hence, straightforward continuation of the perturbation pro- 
cedure breaks down. 


The difficulty arises from the fact that in the first perturbation, 
the 3 / 2 -power branch point in g arises at the basic cone (0 = b) rather 
than at its actual location on the ogive surface. Although this dis- 
crepancy has no effect upon the first perturbation, it is compounded in 
subsequent terms so as to be catastrophic. The remedy is to choose a 
slightly strained coordinate system such that for each term the 3/2-power 
branch point appears precisely at the body surface. Both the difficulty 
and the remedy are just those considered by Lighthill in his discussion 
of a technique for rendering approximate solutions uniformly valid 
(refs. 27 and 28). As in the previous examples, the solution proceeds, 
in effect, by assuming a given shock wave and determining the corre- 
sponding body shape. Therefore, the required straining of coordinates 
is not known at the outset, but must be determined to successively 
higher accuracy as the solution progresses; this is characteristic of 
Lighthill *s technique. 


With this modification, the perturbation procedure can, in principle, 
be continued Indefinitely. It can therefore be concluded that an analytic 
body of revolution is accompanied by an analytic attached bow shock wave 
at supersonic speeds, and an analytic pressure distribution. 


Uniformly valid equations of motion .- Let the body be given in 
phy 8 ical coordinates by 


= T ^bx + |cx^ + ^dx® + . . . 


and the corresponding shock wave by 

r. = T X + izcx® + i(md + nc^)x® + ... 


(67) 


( 68 ) 
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Now introduce a slightly strained radial coordinate r such that the 
body surface is given by r = bx. The simplest choice is 

r = f - - g< 3 x'^ + ... (69) 

The procedure which led to equation ( 57 ) gives for the entropy 
function behind the shock wave 


to(i|f) = 0 )q 



iZCi/2ijf + [u)i(md+nc^) + to^Z^c^] \|f + 



where cJq and 0)i are given by equation ( 40 b)^ and 

(O2 = Ml+i) 

(Ko2_i)[2+(7-i)ko2] 


(70a)' 


(TOb) 


The stream function has the form 


\|f = x^f(e) - cx^(e) - x'^[c^j(0) + dh(e)] + ... (71) 


where ? is the nearly conical coordinate r/x. The differential 
equations which result from substituting this series into the equation 
of motion are simplified by setting 


g('^) = lg(8) - -if '(8^) 

S(8^) = mh(8^) - (S) 

6 


} (72) 


■^(^) = Z^j(e) + nh(8T) + |zg-(^) - if-(8r) J 

(The functions g, h, and j thus introduced axe those which would 
appear in place of g, h, and ^ if no straining of the coordinate system 
had been luidertaken. ) Then the differential equations (and boundary 
conditions) for f and g, as functions of the strained variable 9, are 
(as implied by the common notation) found to be just those solved in 
the previous sections, where f and g were functions of the unstrained 
variable 9. The differential equations for the new functions h and j 
are 


Dh" = E + Fh + Gh' 7 
DJ" = I + Fj + Gj' j 


( 73 a) 
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where D is given hy equation ( 58 b), and 

fi(r+i) 

E = 

0 ( 7 - 1 ) 

F = 4 (f'^ + 4 ff") 


yf[ f- . |1 


+ f 


G = 7^0 


.,(7+1) r 


>(7-1) 


2i£ - (7+1) — 
e f ' 


- 12ff' 


I = 7Uq 


m( 7 +i) 

n(7-i) 


0) 


2 ^,2 


(7+1) 


-s' 



g"- Si 


(Oof* + 


7(7+1) -g» 


2 -I 


(0 




2(7+1) IzilZsl . 1 
f'/f 


i f" \ + f 'g' 

9 


7+3 + f'/f 
7 27 (-g'/g)^ 


3g^f" - 2gg'f - 2(g'2 _ 2gg")f 


(73b) 


Boundary conditions .- Again, ^expressing conditions just behind the 
shock wave in terms of those at 0=1 by Taylor series expansion and 
simplifying with the aid of the differential equations gives 

\ 

h(l) = i g(l) 


h'(l) 


1 

3 


g.(l) , 

2+(7-1)Kq2 J 


> (74a) 


j(i) = - 1 


2g'(l) + 


2+(7-1)Kq^ J 


J'U) = 1 - i g"(l) - lt( 7 *l)K / 


J 
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where f*(l) is given by equation (51b), g(l) and g'(l) by equation (59),' 
and 


2+(7-1)k^2 

^ 2Kq^E27ICo^-(7-1)] [2+(37+1)Kq^] 

(7+1)(Ko^-1) [2+(7-1)Ko^J2 

g"(l) = 


)(74b) 


K 2 ^(7-1 ) (27+3 )+^(273+37^+5 )kq^-3(7-1 ) 7^+(67+1)kq‘^+2 (7-1 ) (57+3 ) kcP 

° (7+1)(ko^-1) [2+(7-1)Kq2]® 


Behavior near body surface .- The fact that f(?T) is analytic 

at ? = b implies that the full solution for h is also analytic, and 
that the solution of the homogeneous equation for J is analytic. The 

coefficient I is proportional to (?f - b)"®'^^, so that a particular 

integral for j is proportional to (? - b)^^^. For purposes of computa- 
tion it is important to separate the regular part of I from the singu- 
lar part, because either may predominate for the closest practical 
approach to the singularity, depending upon the value of the similarity 
parameter. Furthermore, the accuracy of joining the step-by-step solu- 
tion for j to the series is increased by treating not j itself but 
the combination j + g'/2Z. Hence, the series employed is 


2 Z 


a/2Wi m 

b L 8 

7 

Sbgb’ 

1+ € + (7_J_ ) vs^ + . 

fb' 

L V 2n / 

J(b) 

1+2 Ae^ + i2 Ae® + ... 


tOi 


Wn-WoA .._V2 


— • • « 


J'(b)(0-b) 




(75a) 


where 


fb = f(b), etc., A and e are defined by equation (64b), and 
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M = 


V = 


Nmnerical integration 

have been integrated numerically in the manner outlined previously, 
with 7 = 1.405. Because h represents actually (like g) o n ly a first 
Perturbation of the basic conical flower and is furthermore regular near 
the body, it is readily determined with ample accuracy. On the other 
hand, in the ecjuation for j, the coefficient I is so strongly singu- 
lar that it was fornid necessary to replace simple step-by-step integra- 
tion by the more laborious five-term procedure of Milne (ref. 29, 

P* 1^2). The coefficients and boundary values are also considerably 
more difficult to calculate, so the integration of j has been limited 
to three values of the similarity parameter, whereas h has been found 
for four values. 


gb 

2/2 U)2 (bf^ ' 

37 gb 

The differential equations for h and i 


(75b) 


The accuracy of the solution for j suffers from the facts that it 
depends upon the accuracy of the preceding solution for g, that one 
coefficient in the differential equation is strongly singular, and that 
the results of physical interest are found as differences of nearly equal 
quantities. Consequently, although results derived from the fxmctions f, 
g, and h are probably reliable to three or four significant figures, ^ 

those derived from j are perhaps not reliable to more than two. 


Body shape and pressvire distribution .- The parameters m and n, 
which relate the shape of the shock wave to that of the body, are found, 
by requiring the stream function to vanish at the surface, to be given 


. m 



r 


n 



Ob J 


> 


(76) 


The surface pressure coefficient is given by 



T 


2 


2 

7Ko^ 




(77) 


It corresponds to a pointed body with zero initial curvature; 
compare reference 24, part III. ^ 
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Thus it is found that on the surface of a body described in actual 
coordinates by 


r = Rq'x + i Ro"x^ + ^ + ••• 


( 78 ) 


the pressure coefficient is given by 


P _ p 1 2 / 

^TV ~ -RO 




Rr 


+ Ro'Ro" 


L8(Ro’Ro"x) J 


X + r ^"2 r 




21 ° L8 ^(Ro"x)^ 




Ro’Ro 




L8(Ro'x)a(Ro'"x) J 


X + ... 


(79a) 


where 


R 


= 2 


■ f vy L_ 

_ b^ V b / r(MRr,' 


r(MRo')^ J 


8(R. 




^ /^bV /~gb*\ 

^o'Ro"x) b W y V gb y 


2 TWc 


(^vy / 

8(Rq'x)8(Rq'”x) 3 b ^ b \ 






>( 79 b) 


8(Ro"x)^ 


°V »9 W 85= 2b= J5 t^O) 


Numerical values of h(b), h'(b), j(b), and j'(b) are listed in 
the following table, together with the resulting values of the parameters 
m and n which relate the shape of the shock wave to that of the body, 
and the two functions which give the initial curvature of the surface 
pressure distribution. 
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pRo’ 

h(b) 

h’(b) 

J(b) 

J'(b) 

m 

n 



8(Ro’’x)2 

8(Ro'x) 8(Ro’"x) 

0.3765 

.6599 

1.150 

09 

l 4 . 3 i^ 

1.752 

• 93 ‘H 

1.166 

-36.40 

- 5.401 

-3.665 

- 6.948 

-7.34 

1.94 

7.63 

43.7 

5.30 

-27.4 

0.0051 

.0803 

.2867 

.8031 

0.583 

.868 

-.901 

4.04 

5.20 

13.7 

5.694 

4.736 

4.561 

5.332 


The two pressure-curvature functions are plotted in sketches (n) and (o). 
(Curves have been faired through the three calculated points by analogy 
with the results of the cone -expans ion approximation discussed in the 
following section.) 





Sketch (n).- First term in ini- 
tial pressure curvature on 
ogive of revolution. 



Sketch (o).- Second term in ini- 
tial pressure curvature on 
ogive of revolution. 


Fiirther Approximations 


The theory discussed heretofore is the simplest which retains all 
the essential features of hypersonic flow, so that its solutions approach 
exactness as the thickness tends toward zero. Further approximations, 
although desirable for facilitating solution, will introduce errors 
whose nature may be more obsciure. In the case of plane flow, however, 
further approximations exist which are so simple and accurate that the 
problem may be considered solved for practical purposes (cases ( 6 ) and ( 7 ) 
below). These and other approximations will now be considered for three- 
dimensional shapes, in comparison with the solutions already given. 

Most of these approximations are useful outside the limits of hyper- 
sonic small-disturbance theory. However, we shall consider them here 
only as they are reduced to small-disturbance form, so that they actually 
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represent approximations beyond those already made. For example, the 
well-known shock-expansion method will be considered only in its hyper- 
sonic small-disturbance form (ref. 17). 

The following additional approximations will be considered: 

(1) Linearized theory, second-order theory, etc. 

(2) Newtonian impact theory 

(3) Newtonian theory plus centrifugal forces 

(4) 7=1 

(5) Cone-expansion approximation 

(6) Tangent-cone approximation 

(7) Compression-layer approximation 

Linearized theory, etc .- The breakdown of linearized theory serves 
almost as a definition of hypersonic flow. Hence, the most that can be 
expected of linearized theory, second-order theoiy, etc., is that t'hey 
penetrate somewhat into the lower end of the hypersonic range. 

For plane flow Donov (ref. 30, pp. 90-91) has determined the fourth- 
order solution. Reduced to hypersonic small-disturbeince form, his result 
for surface pressure coefficient on a single airfoil may be written as 

= 2K + Z±1 K^+(y+l)(37-?) + + 3-7-37^+/^ ir 4 + 

p 2 6 48 48 32 ^ 


3 :77-77^1-37^ Ko=k + 

96 


(7+l)^(37-^) 

128 


Ko®Ko'x 


+ 


(80) 



Sketch (p).- Further approximations 
to hypersonic small-disturbance 
theory; pressure on cone, 7=1.405 


where K is the local similarity 
parameter (M times local surface 
slope), Ko its value at the lead- 
ing edge, and Kg' its initial rate 
of change. Even in this reduced 
form nothing is known of the range 
of convergence of the series or, 
indeed, whether it converges at all. 
However, for a single wedge the 
solution is known in closed form 
from equation ( 3^) • Hence , it ia 
seen that in this special case the 
series is convergent for 
K s M6 < 4/(7-t-l), which is 1.67 
for air. 

For cones the linearized and 
second-order solutions (ref 31> 
p. 11)^ as reduced to hypersonic 
small-disturhance form, are shown 
in sketch (p). 
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Nevtonian Impact theory »- Assuming that fluid particles lose their 
normal momentum on impact with the surface leads to a prediction of pres- 
sures proportional to the square of the sine of the angle of inclination 
or, in the small-disturbance approximation 

= 2 (81) 


R'^ 


wherever the slope is positive, and zero elsewhere. According to equa- 
tion (3ll-) the actual value for a wedge falls only to 2.1+ at infinite 
Mach number (with 7 = 7/5)^ so that the approximation is poor for plane 
flow. It is more satisfactory for fusiform shapes such as a cone 
(sketch (p)), for which the actual value at M = « (with 7 = 7/5) is 

2.09. 


Newtonian plus centrifugal forces .- Newtonian impact theory has been 
improved by including the centrifugal pressure gradient through the layer 
of fluid streaming over the body (refs. 32 and 33). The result is pre- 
cisely the limit of the full theory as M->- » and 7-> 1. In the small- 
disturbance approximation (ref. 33), it gives for plane flow 


Cp = 2(R'^ + RR") (82) 


and for axially symmetric flow 


Cp = 2 R '2 + RR" (83) 


In both cases it is to be understood 
that negative values are to be 
replaced by zero. Sketch (q) shows 
that the improvement due to includ- 
ing centrifugal effects is appreci- 
able for the initial pressure gradi- 
ent on an ogive of revolution. 

7=1.- It has just been seen 
that on fusiform shapes near 
M = «o, the surface pressure is 
insensitive to the value of 7. 

At the other end of the hypersonic 
range, linearized theory is inde- 
pendent of 7. These two extremes 
suggest that a close approximation 
throughout the hypersonic range 
may be found by setting 7=1 
(and this is particularly tinxe 
since in a real gas 7 



•333 .5 I 00 

Sketch (q).- Further approxi- 
mations to hypersonic small- 
disturbance theory; Initial 
pressure gradient on ogive 
of revolution, 7 = 1.405. 
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approaches 1 at high temperatures). This choice simplifies the theory 
by rendering it effectively isentropic; that is, although shock waves 
produce entropy jumps, entropy does not appear in the pressure-density 
relation and is therefore absent from the problem. 

This approximation has been tested by computing the hypersonic small- 
disturbance solution for a cone with 7=1. The results corresponding 
to those tabulated on page 27 are shown in the following table (including 
the known value at M = 00 ) . 


Mt 

5 /t 

p& 

f'(b)/b 


1.04 

0.3912 

0.4069 

1.255 

3.076 

1.15 

.5609 

.6450 

1.546 

2.627 

1.50 

.7647 

1.147 

2.498 

2.277 

00 

1 

00 

00 

2 


Sketch (p) shows close agreement with the results for 7 = 1.405^ the 
discrepancy being, indeed, less than that due to the thickness of a 10° 
semi vertex angle (cf. sketch (k)). 

"Cone -expans ion" approximation . - The shock-expansion method for 
plane flow, which neglects disturbances reflected from the bow wave, has 
recently been shown to yield good accuracy at all supersonic speeds away 
from the transonic zone (ref. 34). 

A more surprising discovery is that an analogous procedure yields a 
reasonable approximation for certain three-dimensional shapes in hyper- 
sonic flow. In this "cone-expansion" method the flow behind the tip of 
a pointed body is approximated by a Prandtl-Meyer expansion (refs. 35 
and 36). The accuracy of this approximation is indicated by the compar- 
ison shown in sketch (q.) for the initial pressure gradient on an ogive of 
revolution. 

Tangent -cone approximation . - Newtonian Impact theory predicts pres- 
sures depending only upon the local slope. This suggests approximating 
the pressure at each point of a body by that on a locally tangent cone 
or wedge at the same Mach number. For plane flow this gives equation (35)^ 
which yields good accuracy. For bodies of revolution, sketch (q) gives 
an indication of the accuracy obtainable. 

Compression-layer approximation . - In the upper end of the hypersonic 
range the bow shock lies close to the body (if the body slope is posi- 
tive). This suggests making an approximation somewhat analogous to that 
of the Prandtl boundaiy-layer theory, assuming that the layer of com- 
pressed fluid between the body and shock is very thin. 

For example, assume that the shock wave lies so close to a circiilar 
cone that a linear variation is adequate to describe the stream function. 
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Then according to eq,uation 


f( 0 ) 


Requiring this to vanish at the surface gives as the ratio of cone angle 
to shock-wave angle 

^ _ 6 _ (r-^3)Ko^-g 

’’ 2(7+1)I'o* 

This result has been derived by Lees (ref. 37 ). At infinite Mach num- 
ber with 7 = 1.405 it gives O.916 compared with the true value of 0 . 9 l 4 , 
and Lees shows that it is accurate even in the lower end of the hyper- 
sonic range. However, for the corresponding surface pressure coeffi- 
cient, the range of good approximation is much smaller, as shown in 
sketch (p). 


(51) the stream function is given by 


(7-H)Kq^ 

2+(7-1)Kq" 


( 1 - 0 ) 


( 84 ) 


Ames Aeronautical Laboratory 

National Advisory Committee 
Moffett Field, Calif., 


for Aeronautics 
Mar. 18, 1954 


kk 


MCA TN 3173 


A^B^C^D^E^ 

F,G,H,I 

B(x,y,z) 

b(x) 

b,c,d 

Cp 

^ j 
I 

in,n 

M 

P 

R(x) 

r 

S(x,y,z) 

s(x) 

t 

u,v,w 


APPENDIX A 
PRINCIPAL SYMBOLS 

coefficients of differential equations. (See eqs. ( 58 ) 

and (73).) 

fimction defining body shape 

reduced radius (or ordinate) of axially symmetric (or 
plane) body 

coefficients in series expansion for radius (or ordinate) 
of body (See eqs. (36) and (67).) 


pressure coefficient, 2_ 

1- p 

- P«u« 

functions in series expansion for stream function (See 
eqs. (39) and (72).) 

initial ratio of shock-wave curvature to body curvature 

coefficients relating shapes of shock wave and body 
(See eq. (68) . ) 

free-stream Mach number 

pressinre 

radius (or ordinate) of axially symmetric (or plane) 
ogive 

radius (or ordinate) in cylindrical (or plane Cartesian) 
coordinates 

function defining shape of system of shock waves 

radius (or ordinate) of shock wave attached to axially 
symmetric (or plane) body 

time 

velocity components in Cartesian or cylindrical 
coordinates 
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1^5 


x,y,z Cartesian coordinates with x in stresunwise direction 

3 a/ - 1 

7 adiabatic exponent 

6 semi vertex angle of wedge or cone 


(e - b) 

b 


e 


conical variable, 


X 


k(x) 


A 


aiaxi'liary hypersonic similarity parameter based upon local 
slope of shock wave, Mts'(x) 




-1 




T^(-Sh*) 

Sb 


2»/2(az ) 

37 gb 


P density 

o constant which is zero for plane flow, unity for axially 

symmetric flow 

T thickness parameter of body; in examples, initial slope of 

shock wave 


^ stream function for plane or axially symmetric flow 

(See eq. ( 23 ).) 

w entropy function (See eq. (25).) 

(0o»W3_,c>)2 terms in series expansion for u (See eqs. (4o) and ( 70 ).) 

[ ] denote jump in quantity through shock wave 
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kS 


( ) 

/v» 

( ) 

( )' 

( )o 
( >b 

( h 
( ). 


reduced form (See eqs. (8).) 

form associated with strained coordinates. (See eqs. (69) 
and (71.)) 

derivative with respect to argument 

value at tip of pointed body- 

value at surface of body (or at 6 ,G = b) 

value at shock wave 

value in free stream 


APPENDIX B 


CONNECTION BETWEEN HYPERSONIC AND LINEARIZED 
SUPERSONIC SIMILITUDE 


The similitude for linearized supersonic flow is now well understood, 
having been first correctly stated by Gothert in reference I5 (for the 
analogous case of linearized subsonic flow) . This similarity rule implies 
that the reduced coordinates x,y, and z of equation (8a) may again be 
introduced, and that then the reduced flow quantities. 







(Bla) 


(Bib) 


depend only on^the reduced coordinates and the supersonic similarity 
prameter px. The error- in the theory and associated similarity rule 
is 0(t/)3), in general. It may be emphasized that here, as in all the 
similarity rules, the choice of reduced yariables is by no means unique: 
an unlimited number of equivalent forms can be produced, for example, by 
multiplying each reduced variable by, or adding to it, any constant 
multiple of powers of the similarity parameter. The particular forms 
adopted here were chosen to correspond as closely as possible to their 
hypersonic counterparts in equation (8b), so as to facilitate the follow- 
ing argument . 

For Mach numbers so large that M is effectively equal to p, 
these results agree with those for the hypersonic case,i3 and this'was 
pointed out by Tsien (ref. 5). However, it is more fruitful to reverse 


_ _Here, in contrast with the hypersonic case, the reduced variables 

u and p are definitely not 0 (l) as T — >0 (for fixed M) . 

^^With p" = 1 + ^ and p = ^ + 
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the argument, and observe that the hypersonic similitude, Just as it 
stands, is entirely consistent with the linearized similitude. This is 
immediately apparent for the reduced velocities u, v, and w which (as 
implied by the common notation) have identical forms in the two cases. 
They differ only in depending upon different parameters, but in hyper- 
sonic flow M and p are interchangeable to within the accuracy 0 (t^) 
of the theory, so that Mr can be replaced by Pt to complete the 
correspondence. For the pressure and density, the hypersonic theory 
(eq,. (8b)) shows that 


1 

_ 1 


7M^T^ 


/ 7M^T^ 

gf 

VP« 

) = 


(B2a) 


Again utilizing the fact that M and p are interchangeable in the 
hypersonic range, these can be rewritten 


— ~ ^ 2^ “ P 




(^0 


(B 2 b) 


= P - 1 


= P 


where the final forms depend upon the parameter pT and are therefore 
(as implied by the notation) identical with their linearized counter- 
parts in equation (21b). Thus the correspondence is complete. 


This means that the hypersonic small-disturbance theory, when 
properly Interpreted according to. the linearized supersonic similitude, 
yields a first-order solution at all speeds above the transonic zone. 
The reduced problem of equations (lO) to (13) is solved for a given 
value of the parameter Mr, and then with Mr replaced by Pt is 
interpreted in terms of physical variables according to 

u = u„ [1 + T^(x,y,z; pr)] 1 


V = U„TV 


w = u„ TW 


p = [1 . (p - ] = p„(rM^r2p . 


> (B3) 


P = Po 


1 + i!, (p- 1) 


= Po 
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